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Abstract 

A theoretical investigation has been carried out to explore the use of electrically conducting fluids as lubricants in circular 
and elliptical hydrostatic thrust pad bearings subjected to transverse magnetic field. Flow of electrically conducting fluid in 
the presence of external transverse magnetic field is represented by the addition of Lorentz force in Navier-Stokes 
equation. A MATLAB source code based on finite element formulation of the modified Reynolds equation is developed 
for obtaining static and dynamic characteristics parameters of the bearing. Influence of aspect ratio (major/minor-axis) 
and magnetic field (Hartmann number) are analyzed on the performance of elliptical pad bearing. Application of elec- 
trically conducting lubricant in the elliptical pad is observed to provide higher values of fluid film stiffness (44.93%) and 


damping (139.37%) coefficients as compared to the circular pad operating with Newtonian lubricant. 


Keywords 


Hydrostatic thrust pad bearing, electrically conducting lubricant, magnetic field, elliptical pad, finite element method 


Date received: 16 August 2017; accepted: 13 December 2017 


Introduction 


Hydrostatic thrust pad bearings are used for support- 
ing the axial load in the rotating machinery and struc- 
tures. These bearing are employed in the applications 
ranging from small precision machine tools to sup- 
porting some of the heaviest machinery and heavy 
structures such as turbomachinery, vertical pumps, 
telescopes, observatory domes, etc. Use of the hydro- 
static thrust pad bearings enables complete separation 
of matting surfaces during start-up and stoppage of 
machinery. If needed, a complete separation of bear- 
ing surfaces by thick fluid film is even possible at 
dead/zero speed of runner. These classes of bearings 
offers high fluid film stiffness and damping coefficients 
as well as exact positioning and control in engineering 
applications vis-a-vis to the hydrodynamic bearings. 
Owing to wider applications domain, tremendous 
research activities have been carried out in the 
hydrostatic thrust pad bearing in the last few decades. 
In the beginning, Fuller’ investigated the squeeze film 
performance of the hydrostatic bearings. Later, 
Malanoski and Loeb? investigated the influence of 
restricting devices i.e. capillary, orifice, and constant 
flow valve on the fluid film stiffness of the hydrostatic 
thrust pad bearing. Theoretical? investigation was 
reported concerning the influence of elastic distortion 
of the bearing liner on the performance of the 


hydrostatic thrust pad bearings. The elastohydrostatic 
behavior of the bearing liner in heavily loaded appli- 
cations was observed to have beneficial effect on the 
bearing performance. Osman et al.* used a specially 
designed test rig to assess the performance of circular 
and sectorial recessed hydrostatic thrust pad bearing. 
The authors concluded that the recess shape, size, and 
location significantly alters the bearing performance. 
Maher et al.* presented an alternate approach to 
predict the performance of the noncircular recessed 
hydrostatic thrust pad bearings. Pressure distribution 
in the rectangular? and elliptical® recess hydrostatic 
thrust pad bearing were obtained by employing 
Riemann’s mapping theorem and Jacobian elliptic 
functions to map these recess shapes to circular 
recess in the hydrostatic thrust pad bearings. 
Lubricants used in the aforementioned studies!” 
were assumed to be classical Newtonian fluid i.e. 
having linear relationship between shear stress and 
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strain rate. However, lubricants used in typical indus- 
trial applications are generally blended with 
certain polymeric additive packages to improve the 
tribological performance of the bearing systems. 
These additives are also desirable to limit the vari- 
ation in lubricant viscosity at elevated operating 
conditions. Researchers in the past have theoretically 
and experimentally investigated the squeeze film oper- 
ation of the thrust bearing lubricated with non- 
Newtonian fluid such as power law fluid,’ couple 
stress fluid, ° '? and Carreau model fluid!? in various 
classes of fluid film bearings. Recently, Yadav and 
Sharma? theoretically investigated the performance 
of hybrid thrust pad bearing with micro-patterns on 
the bearing surfaces, operating with non-Newtonian 
power law lubricants. A significant reduction in the 
power losses was observed by using micro-textured 
surfaces and non-Newtonian power law lubricants in 
the hybrid thrust pad bearings. All abovementioned 
studies reported significant improvement in the bear- 
ing performance owing to the use of non-Newtonian 
lubricants. One particular non-Newtonian lubricant 
that has received considerable attention over the last 
few decades is the electrically conducting fluid. Use of 
electrically conducting fluid as lubricant in the pres- 
ence of external magnetic field stabilizes the lubricant 
viscosity variation with the temperature. A number of 
researchers'**° have focused their attention on the 
application of this lubricant in different class of bear- 
ings. Application of liquid metals in hydrostatic thrust 
pad and slider bearings was investigated by Elco 
and Hughes.” Further, squeeze film characteristics 
of parallel circular and rectangular plates lubricated 
with electrically conducting fluid in the presence 
of magnetic field was investigated by Kuzma.'° The 
abovementioned theoretical studies'* '® predicted an 
increase in the load-carrying capacity and time of 
squeeze film approach, due to the use of electrically 
conducting lubricants (ECL). Singh and Malik!’ 
derived two-dimensional modified Reynolds equation 
for magnetohydrodynamic (MHD) journal bearings. 
Solution of this modified MHD Reynolds equation 
was presented in terms of double series. Further, 
study from Shukla and Isa!’ predicted an increase in 
the load-carrying capacity of externally pressurized 
porous thrust bearing on increasing the permeability 
and strength of the magnetic field. Experimental stu- 
dies were carried out exploring the use of Gallium 
(Ga) based liquid metal alloys as the lubricant in 
fluid film journal bearing by Gerkema.'? Influence 
of liquid metal lubricant properties was studied 
on the loci of shaft center of plain and spiral groove 
journal bearing. Gerkema proposed the use of 
Ga-based alloy in the spiral groove journal bearing 
in the nonoxidizing condition. The lubricants with 
higher electrical conductivity i.e. ECLs and liquid 
metal alloys offer better rheological performance vis- 
a-vis conventional nonconducting lubricant in the 
fluid film bearings operating at elevated temperatures. 


ECLs are quite important for the bearings operating 
in the nuclear and space engineering applica- 
tion.'”!?*! Liquid Ga-based metal alloys have the 
ability to wet most metallic and non-metallic surfaces 
and offer excellent lubricity in high temperature bear- 
ing applications.” Researchers also explored 
the use of ionic liquid (IL) lubricant in the elevated 
temperature application. High electrical and thermal 
conductivity of ILs ensures the thermal stability and 
desirable lubricant properties under such conditions. 
MobilTM DTE 932 GT is an example of commer- 
cially available high performance turbine oil? 
having good electrical conductivity. Besides this, the 
use of magnetic field has great importance in engin- 
eering applications such as MHD power generation, 
MHD accelerator, purification of crude oil, MHD 
propulsion, etc. Theoretical studies*® have been per- 
formed on the slider bearings operating with couple 
stress fluid in the presence of magnetic field to obtain 
optimum load-carrying capacity and input/output 
film height. Influence of surface roughness on the 
MHD squeeze performance of porous parallel rect- 
angular plate was investigated by Kesavan et al.” 
It has been observed that the presence of magnetic 
field and negatively skewed surface roughness is 
beneficial in enhancing the load-carrying capacity of 
parallel rectangular plates. Influence of the geometric 
shape of plates on the MHD squeeze film behavior 
of finite parallel plate lubricated with conducting 
couple stress fluids was analyzed by Fatima et al.78 
The authors concluded that the use of circular parallel 
plates yields highest transient load-carrying capacity. 
Bompos et al.” investigated the use of magnetorheo- 
logical (MRF) and nano-magnetorheological fluids 
(NMRF) in the rotor dynamic behavior of journal 
bearings. Result obtained from CFD simulation in 
the study proposed the use of NMRF for improving 
the static and dynamic performance of journal bear- 
ings. A recent study from Kumar and Sharma?’ 
reported a substantial enhancement in the dynamic 
characteristics of hydrostatic thrust pad bearing lubri- 
cated with MHD-couple stress fluids. Lubricant used 
in the common industrial applications possesses very 
high electrical resistance. Prasad investigated the per- 
formance of the pivoted pad thrust bearing?!“ and 
rolling element bearings?!” operating in the electrical 
environment. The author reported that the presence 
of electric field accelerates the wear and lubricant 
deterioration, which will subsequently reduce the life 
of bearings. This is due to the discharge of accumu- 
lated charge between bearing surfaces, in the form of 
sparks resulting/initiating surface fatigue. This prob- 
lem can be overcome by the use of lubricant having 
high electrical conductivity and employing insulated 
bearing surfaces. Under such operating conditions, 
the presence of external electric/magnetic field does 
not lead to the accumulation of charge over the bear- 
ings surface, due to high electrical conductivity of the 
lubricant. 
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A thorough review of the available literature 
indicates that studies associated with ECL are 
mostly confined either to hydrodynamic lubrication 
regime or uncompensated hydrostatic parallel plates. 
Coupling of compensator/restrictor flow equation 
with Reynolds equation is essential to accurately pre- 
dict the performance of thrust pad bearings, which is 
missing in majority of these studies. It may be realized 
that the restrictor is a very important element of a 
hydrostatic bearing system for adjusting the variation 
in the external load coming to the bearing corres- 
ponding to the variation in the fluid film thickness. 
A study of the hydrostatic bearing without a restrictor 
is impractical. In addition to the restrictor, the per- 
formance of a hydrostatic thrust pad bearing can be 
effectively controlled by appropriate selection of geo- 
metric shape of thrust pads. In the past, majority of 
research works on these bearings has been focused 
on obtaining the optimum shape of recess, whereas 
the influence of the outer shape of pad on the bearing 
performance is relatively unexamined. Therefore, this 
study is planned to explore the use of ECL and ellip- 
tical shape pads in the hydrostatic thrust pad bear- 
ings. Thus, th objective of this study is to investigate 
the performance of a capillary compensated elliptical 
hydrostatic thrust pad bearing operating with the 
ECL in the presence of the magnetic field. A paramet- 
ric study is also carried out to examine the influence of 
aspect ratio (major to minor axis) of the elliptical pad 
and strength of magnetic field (Hartmann number) on 
the bearing performance. 


Mathematical formulation 


Figure 1 depicts the general layout of the hydrostatic 
lubricating system under investigation. The flow in the 
bearing is assumed to be laminar, incompressible, iso- 
thermal, and isoviscous. The incompressible electric- 
ally conducting fluid is supplied to the recess of the 
bearing pad via capillary restrictor in the presence of 
uniform external magnetic field (B,) along the fluid 
film thickness direction. Using the conditions of 
incompressible fluids and steady-state flow: # — 0, 
V.¥ =0 and neglecting the inertial forces, body 
forces (except Lorentz force), couple stress, and 
body couples, the Navier-Stokes and flow continuity 


equations are expressed as 


O=-Vpt p+ J x B (1) 


> ou v ðw 


=0 (2) 


> f 
Here, F; represents the Lorentz force experienced 
by a charge particle in the presence of the magnetic 
field (B). Lorentz force is a body force and is 
> > > =>. 
expressed as: F; = J x B. Here, J is the current 
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Figure |. Schematic diagram of an elliptical hydrostatic thrust 
pad bearing system. 


density and defined by: J= Ae +9 x B). 
Momentum equation (1) in simplified form along x, 
z, and y directions can be written as 


p > 3 
u P — o Bu = a +ocE-B, (3) 
ew 2 op 
Page One a, ee (4) 
op 
—=0 3 
© 


E, and E, are the components of induced electric 
field in the x and z directions respectively. Assuming 
the bearing surfaces to be the perfect insulators, the 
electric field in equations (3) and (4) can be obtained 
by using the condition of zero net current flow i.e. 


h 
f (E, + uB,)dy = 0 (6) 
0 


h 
(Ex — wBo)dy = 0 (7) 
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Boundary conditions used for obtaining velocity 
profile are expressed as 


oh 


Ono slip’ S ‘Ot (8) 


Aty=0; (w= U; w 


Aty=h; (u=0; v=0; w=ONostip (9) 


Integrating equations (3) and (4) using the above- 
mentioned boundary condition (equations (8) and 
(9)), and simultaneously solving equations (6) and 
(7) for E, and E., the final expression for the velocity 
of fluid in x and z directions is expressed as 


U (sinn(42) = sinh(“4—) ) 
"=z ( sinh (#4) 
hoh op (sinn (4) sinh(42) sinh (44-2) ) 


2H dx (cosh (44) = 1) 


(10) 


ny (en) — (2) = sE) 


2uH dx (cosn(44) — 1) 


wW = 


(11) 


Here, H is the Hartmann number denoting the rela- 
tive strength of the electromagnetic force to the viscous 
force, expressed as: H = B,h,(o/1)'/* where o denotes 
the electrical conductivity of the lubricant. Substituting 
the value of velocities from equations (10) and (11) 
into the continuity equation i.e. equation (2) 


ge (oH) +5 a(r oh, mn) 
Ox 0. 


oh oh da 
1 h 
= UU + 12 
where 
6h? Hh Hh 
o(h, H) = Ae corn (2) 2) (13) 


is the MHD function 


Equation (12) is the modified Reynolds equation 
for ECL in the presence of external magnetic field 
(B,), and hydrodynamic terms in the above equation 
can be ignored to study the performance of the hydro- 
static thrust pad bearings. In the limiting cases, the 
MHD function (equation (13)) reduces to unity 
dim (h, H) = 1) i.e. the Newtonian lubricant (NL). 


Using the following nondimensional parameters 


7 sa uf 
a=% > Ah ; ; gee) 
ro ro lo Ps 
14 
- Ips Uro (14) 


Dimensionalless form of the modified Reynolds 
equation is expressed as 


op op əh oh 
da (F gh, H) 2) (F gh, H) 7) -= VaT 


(15) 


where 


oh, H) = (< : ) (ancon (> ZA -2) (16) 
? FP 


Finite element method (FEM) approach has been 
used to solve the two-dimensional modified Reynolds 
equation (equation (15)). Thrust pad surface is discre- 
tized using the four-noded bilinear isotropic quadri- 
lateral element. The fluid film pressure within the 
four-noded quadrilateral elements is approximated by 


4 
p= XO iN; (17) 
fl 


Here, N; is the elemental nodal shape function for 
the four-noded quadrilateral element. Using the 
approximated value of pressure (equation (17)), the 
residue of the Reynolds equation is obtained as 


a f7 i,m a 
Ja h g(h, aa (2? j 


afa —-..af<~ ðh əh 
ral? al.) — 22° — = = Res 
(18) 


Equation (18) presents the weak formulation of 
modified Reynolds equation, with Langrangian inter- 
polation functions (N,) as weight functions 


| bs (N; * Res)dadB = 0 (19) 


After FEM formulation and application of 
Galerkin’s orthogonality conditions on equation 
(19), the elemental system of equation for obtaining 
nodal pressure vector is expressed as 


[Fs] (°} = {Os} + [Re] (20a) 


(o°} = [E] [tO + AR] (200) 


Matrix terms in equation (20a) for eth element in 
the discretized parallel thrust pad surface is expressed 
as follows 


me 3 = ƏN; ON; 3 = ƏN; ON; 
B= ff E Oh, Ha ath g(h, H) J5 aaa 
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Fluidity matrix for eth element 


a- (vam 


Flow matrix for eth element 


D) + (Pan H) z) m | Ndr" 


Ri, = | N,dadB 


Squeeze matrix for eth element 

Elemental matrix equation (equations (20a) and 
(20b)) is generated for all the elements and assembled 
to form global system of equation, expressed below in 
the matrix form 


Fi Fi . F Fin Pi 
Fy Fr Fy Fon p2 
F, Fi : Fj Fin Di 
F, Fi te Fj te ae Pn (21) 
Qı Rt 
Q2 Rh 
=| i |+al — 
i ti 
On, Rt, 


The squeeze term i.e. Ri, in equation (21) can be set 
to zero to obtain the steady-state (0h/dt = 0) pressure 
vector. Equation (22) expressed below incorporate the 
boundary conditions i.e. all the nodes lying on the 
recess boundary are assigned equal pocket pressure 
and summation of the lubricant flow through these 
nodes equals the flow of lubricant from the restrictor. 
This will ensure that the lubricant input flow rate is 
equal to the lubricant flow rate through restrictor or 
lubricant leakage rate at the outer boundary of thrust 
pad. All the nodes lying on outer radius of bearing 
pad are assigned zero pressure value and follows 
Reynolds boundary condition. 


Fy Fy Fj Fin |T pi Qı 
Fy, Fo Fy; Fon P2 Q2 
F; Fi Fj Fin Poi 7 Or 
F, Fo eae Fj vex Fy Pn Qn 
(22) 


The flow through bearing and capillary restrictor is 
expressed as 


QR = CS2(1 — Poc)= Q (23) 

Steady-state load-carrying capacity of the hydro- 
static thrust pad bearings has contribution from the 
pocket area (Poc) and land area due to pressure gradi- 
ents in x and z directions. Thus, the load-carrying 
capacity of bearing is expressed as 


F, = F,| land ara Po leeks area 


j >) ST (Ean ) ain] + dc 


(24) 


Fluid film stiffness and damping coefficient of com- 
pensated hydrostatic thrust pad bearings is obtained 
for dynamic loaded conditions (ah/d¢ 40) and is 
expressed as the partial derivative of the fluid film 
reaction w.r.t. fluid film thickness (h) and squeeze vel- 


ocity (h) respectively. 


oe 2) fle pea Nj) |J|d&dn} + Ace Poe 
ah = oh 
(25) 
- OF 5 z = Poe 
ĉ=z” [ [ 2 vain] E 
oh | ‘=I ðh oh 
(26) 


As stiffness and damping coefficient are dependent 
on the derivatives of pocket pressure and nodal pres- 
sure vector, equation (20) is used to obtain these 
derivatives of the nodal pressure vector. 


a{p} 2 [ET Ee z Rd {p} sy (27) 


Solution procedure 


A source code based on the Gauss-Siedel iterative 
numerical technique is developed to solve the modi- 
fied Reynolds equation (equation (15)) coupled with 
restrictor flow equation (equation (23)). FEA soft- 
ware ANSYS is used to discretize the bearing 
pad surface using four-noded bilinear isotropic quad- 
rilateral element. One circular (O) and four elliptical 
(Ell I-IV) shape thrust pad configurations have been 
investigated in the current study. All thrust pad con- 
figurations have the same pocket and land area. 
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Table I. Geometric and FEM mesh details of thrust pads. 


Circular Elliptical 

Geometric parameters (0) (EII-I) (EIl-II) (EIL-III) (EILIV) 

a | Lt 1.25 1.4285 1.6666 

b l 0.9 0.8 0.7 0.6 

r 0.5 0.5 0.5 0.5 0.5 

a, = a/b | 1.2345 1.5625 2.0407 2.7776 

ma Area of pad 3.141 (r2) 3.141 (rab) 3.141 (xab) 3.141 (rab) 3.141 (rab) 

Nodes 348 526 664 928 1155 

Elements 319 464 588 827 1015 
The geometric and mesh details of different thrust 
pads are shown in Table 1. Moreover, dimensional Table 2. Bearing geometry and operating parameters.?* 
geometric and operating parameters of the bearing : ; 
system are presented in Table 1. The rheological P Dimensional 

a arameters Symbol value 

property of liquid Ga-based alloy (Ga—Jn—Sn), an 
example of commercially available ECL,” have Circular pad outer to 178mm 
been stated in Table 2. The data presented in this radius 
table can be used to obtain the dimensional value of Pocket radius ri 89 mm 
performance parameters of the bearing. To check the Area of bearing (all Ap 99487.76 mm? 
repeatability of the numerical results, a mesh sensitiv- shapes) 
ity test has been performed on the maximum pressure Area of circular pocket Asc 24871.94 mm? 
to achieve an optimum mesh size for the circular Nominal fluid film ho 0.0508 mm 
and elliptical thrust pads. A steep drop in the max- thickness 
imum pressure has been observed on continuous Lubricant supply Ds 20.6843 MPa 
refinement of mesh size. The optimum mesh size pressure 
is achieved when the percentage change in the max- Magnetic field (range) Bo 0-9 Wb/m2 
imum pressure drops below 1%, in the successive Hartmann number H 0-10 
mesh size. Increasing the mesh size above the opti- (range) 
mum value exponentially increases the computational Electrically conducting Liquid gallium(Ga) 51E Indalloy® 
time with any significant gain in the results’ accuracy. lubricant — Indium (In) Number 
The mesh details (Table 1) for the respective thrust — Tin (Sn) alloy 
pads have been obtained following the mesh sensitiv- Composition % Ga — 66.5; 
ity/ grid independence test. In addition to the In — 20.5; 
mesh sensitivity test, Newton—Raphson method is Sn — 13 
employed in the source code to achieve a high con- Lubricant dynamic LL 7.247 mPa-s 
vergence of the order of le~°° for the nodal pressure viscosity (at 20°C) 
vector. The stepwise solution methodology adopted Density p 6499.24 kg/m? 
for computing the steady-state and dynamic perform- Thermal conductivity o 3.46e°° mho/m 


ance parameters (expressed in the preceding section) 
are as follows: 


Step 1: Discretization of the bearing pad surface. 

Step 2: Initialization of parameters (ho, u, p, tolp, H, 
etc.). 

Step 3: Generation of Gauss points for elements. 
Compute elemental matrix (equation (20a)). 

Step 4: Assembly of elemental matrix to form global 
matrix. 

Step 5: Application of boundary condition in global 
matrix equation. 

Step 6: Solving equation (22) for nodal pressures 
using Gauss—Siedel iterative numerical technique. 


Step 7: If the convergence criteria for pressure 
(equation (29)) is satisfied go to step 3, else go to 
step 8. 


Tolp = le“ < max [Apo], (29) 
[Po], 


Step 8: Evaluate static and dynamic performance 
characteristics (equations (24) to (26)). 
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Validation of numerically simulated 
results 


After achieving the mesh sensitivity on pre-set conver- 
gence criterion, the developed source code is validated 
by comparing the numerical result computed from the 
present investigation with the available published 
results.*°* The Reynolds equation, lubricant model, 
geometric, and operating conditions are imported 
from these available experimental and theoretical stu- 
dies for the purpose of validating the program. 
Osman et al.* experimentally and theoretically inves- 
tigated the influence of bearing pad surface roughness 
on the pressure and fluid film thickness of capillary 
compensated hydrostatic circular step thrust pad 
bearing. Bearing pad surface is characterized by flat, 
convex, and concave surface equation. Authors have 
used finite difference method for solving the Reynolds 
equation and obtaining the steady-state performance 
parameters. Using the same geometric and operating 
conditions,* the numerical results are computed 
for flat and convex bearing pad surfaces in accordance 
with the solution methodology described in the 
previous section of the present study. A typical com- 
parison in terms of pressure (Figure 2(a)) and 
film thickness (Figure 2(b)) indicates that the pressure 
distribution and fluid film thickness trends from 
the present study closely follows the curves from 
Osman et al.* A slight variation (<6.5%) between 
the results can be attributed to the different conver- 
gence criteria and mesh size used in both the studies. 
Information regarding the convergence criterion and 
mesh size is not mentioned in the numerical result 
from Osman et al.* In addition to the above, 
the developed program and solution methodology is 
further validated with the theoretical studies con- 
ducted by Fatima et al. presented in Table 3. 
The source code developed in the present study is 
aimed to evaluate the performance of the recessed 
and compensated hydrostatic thrust pad bearing, 
whereas Fatima et al.” evaluated the load-carrying 


capacity for squeeze action between parallel circular 
plates, which are uncompensated and unrecessed. 
The reason for a small variation between the results 
from two studies can be attributed to few assump- 
tions/ modifications made in both studies. First, 
the application of magnetic field on ECLs induces 
an electric field in the fluid film domain. The external 
magnetic field along with this induced electric field 
contributes to Lorentz force acting on fluid particles. 
The reference study** has neglected the induced elec- 
tric field in their analysis. In the present study, influ- 
ence of the induced electric field has been considered 
so as to have a more accurate and realistic analysis. 
This could be the possible reason for the small devi- 
ation between the two results. Secondly, authors 
in the reference study” have not included in their 
analysis any provision (i.e. recess/pocket and restric- 
tor) to facilitate the lubricant supply between the par- 
allel plates. Therefore, to simulate the geometric and 
operating conditions (i.e. nonrecessed and uncompen- 
sated) of the reference study,”® the size of the circular 
recess in the thrust pad for validation purpose is 
reduced in such a way that only five nodes will be 
falling in the recess zone with the value of restrictor 
design parameter being set to a low numerical value 
of 0.01. Thus, the two bearings becomes almost simi- 
lar w.r.t. to geometric and operating conditions. 
The above approximation is attributed to a little vari- 
ation between the results from two studies. Lastly, in 
the current study finite element method is used to 
obtain the performance characteristics of the hydro- 
static thrust pad bearings. Reference study” derived 
closed-form expression to obtain squeeze film charac- 
teristics of parallel plates. Owing to the use of two 
different computational techniques, a minor difference 
between the two results may occur. Aforementioned 
reasons have collectively resulted in the maximum 
deviation of 4.23% between the results from two stu- 
dies, which indicate a reasonably good accordance 
and justify the adequacy of the methodology adopted 
in the present study. 
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Figure 2. (a) Fluid film pressure vs. bearing radius; (b) fluid film thickness vs. load capacity. 
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Table 3. Comparison of dimensionless load-carrying capacity (F,) for different values Hartmann number (H) and 
couple stress parameter (I). 


(I= 0) (i= 0.1) (i = 0.2) (i = 0.3) 

Fatima Fatima Fatima Fatima 
H et al.28 Present et al.7° Present et al.2® Present et al® Present 
0 7.964 | 7.7996 22.7523 22.1655 37.5257 36.249 | 62.0957 59.8717 
l 17.843 17.4948 22.9146 22.5101 37.6889 36.8302 62.2591 60.7376 
2 18.3202 18.7207 23.4007 23.544 38.1773 38.4628 62.7481 63.2117 
3 19.1149 19.5352 24.2108 25.1765 38.9912 40.0915 63.5630 66.2526 


— Circular +++ Elliptical I 


--- Elliptical I = & - Elliptical III 
—< Elliptical IV 


Figure 3. Nondimensional pocket pressure vs. Hartmann number. 


Results and discussions 


The result presented in this section discusses the influ- 
ence of electrically conducting fluids/magnetic field 
(H) and aspect ratio of ellipse (a) on the performance 
of capillary compensated elliptical and circular hydro- 
static thrust pad bearings. Results are presented in 
terms of nondimensional parameter namely pocket 
pressure, lubricant flow rate, load-carrying capacity, 
fluid film stiffness, and damping coefficients. A com- 
parative analysis of the circular and elliptical hydro- 
static thrust pad bearings operating with ECL and NL 
have been presented in the following sub-sections. 


Non-dimensional pocket pressure (poc) 


In this section, the influence of electrically conducting 
fluids/ magnetic field has been discussed on fluid film 
pocket pressure and pressure profile in the land areas 
of circular and elliptical hydrostatic thrust pad bear- 
ings. It can be seen from the results depicted in 
Figures 3 and 4 (A&B), that pocket pressure (Poc) 
and pressure gradients increases with the strength of 


the magnetic field (H) acting on ECL. This phenom- 
enon can be attributed to the presence of Lorentz 
force in the fluid film. The Lorentz force shrinks the 
velocity profile of the fluid between the bearing sur- 
faces. As a consequence, higher quantity of the lubri- 
cant will be accumulated in-between the bearing 
surfaces and this will eventually lead to the rise in 
the pocket pressure and pressure gradients. These 
trends in Poc and pressure gradients are similar for 
both circular (O) and elliptical U-/V) shape hydro- 
static thrust pad bearings. Further, the rate of 
increase in Poe W.r.t. H is comparatively steeper at 
lower values of Hartmann number (H <6). For 
higher values of Hartmann number (H > 10), these 
curves become asymptotic to x-axis. Furthermore, 
the use of circular pad in hydrostatic thrust pad bear- 
ings yields higher Poc compared to the elliptical pad 
bearings. As the aspect ratio of elliptical pad (a,) 
increases, the Poc value tends to decrease. 
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Figure 4. (A) 3D pressure profile for circular (O) and elliptical pad (EI-I!); (B) 2D pressure profile for various pad configuration and 


operating conditions. 
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Contrary to this, the use of ECL in the elliptical 
pad hydrostatic thrust pad bearings produces a higher 
maximum percentage increase in Poc compared to 
the circular pad bearings. The percentage increase in 
the bearing performance parameters (Pees F,, Q, 
S Č) has been evaluated for bearings operating 
with ECL over the NL and is expressed as 


(Rol Fo / S fC )ii=io 
n —(Boc/ F/Q / S / sp 


% Doc = —— = 100 
á MATS 


Circular/Elliptical 


(30) 


The trends and values of the maximum percentage 
increase in Poc for circular and elliptical hydrostatic 
thrust pad bearings are summarized as follows 


13.44% | guy > 9-67%| enm > 8-16%| gnu 
> 7.51% |en_7 > 7.35% |Circutar 


Nondimensional lubricant flow rate (Q) 


The variation of the nondimensional lubricant flow 
rate (Q) for various pad configurations in the hydro- 
static thrust pad bearings is shown in Figure 5. 
The presence of the magnetic field (H) in ECL sig- 
nificantly reduces the lubricant flow requirement. 
Application of the transverse magnetic field on elec- 
trically conducting fluid, flowing through parallel 
plates introduces Lorentz force on the fluid particles. 
This result in a contraction in the velocity profile of 
the fluid and subsequently lubricant flow required 
to support a given load reduces. Based on the numer- 
ical results, the following pattern is observed in the 


lubricant flow_ rate llera < O| en- < Olen 


= PH py m S Ola IV ]a=0—10 : es 
As indicated in Figure 5, a substantial reduction in 
the lubricant flow rate has been noticed owing to the 
use of ECL. Quantitatively, a maximum reduction of 
—48.46% and —48.39% in Ọ have been reported in 
the bearing employing circular (O) and elliptical (J) 
thrust pad due to the use of ECL vis-a-vis to the 
NL. Furthermore, an increase in the aspect ratio of 
elliptical pads has an adverse effect on the lubricant 
flow rate. Reduction in the lubricant flow rate (Q) 
drops from —48.39% to —45.54% as the aspect 
ratio of elliptical pad (a@,) increases from 1.2345 to 
2.7776 (i.e. shifting from pad Ell-I to pad EII-IV). 


Nondimensional load-carrying capacity (F,) 


Figure 6 depicts the dimensionless load-carrying cap- 
acity (F,) as a function of Hartmann number (H). 
The numerical results indicate a substantial enhance- 
ment in the value of F, as the strength of the magnetic 
field/Hartmann number (H) acting on the fluid 
increases. The load-carrying capacity of bearings to 
a large extent depends upon the pressure profile and 
pressure gradients, therefore trends observed in 
the load-carrying capacity (Figure 6) are in line with 
the trends of the pocket pressure (Figure 3) and pres- 
sure gradients (Figure 4). Bearing with circular pads is 
noticed to have higher value load-carrying capacity 
than the bearing with elliptical pads. This phenom- 
enon is observed for all the values of Hartmann 
number (H) i.e. for both NL and ECL. The pressure 
gradients in the land area are strongly influenced by 
the geometric shape of the pad. Therefore, the 
response of the magnetic field on the load-carrying 
capacity of bearings is observed to be actively depend- 
ent on the selection of geometric shape of thrust pad. 
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Figure 5. Nondimensional lubricant flow rate vs. Hartmann number. 
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Figure 6. Nondimensional load-carrying capacity vs. Hartmann number. 
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Figure 7. Nondimensional fluid film stiffness coefficient vs. Hartmann number. 


Although the load-carrying capacity of all thrust pads 
(O and Ell I-IV) used in the hydrostatic thrust bearing 
shows an enhancement by continuously increasing the 
strength of the magnetic field. Yet, the increment in 
the value of F, for various pads configuration is dif- 
ferent and is a function of the aspect ratio/ geometric 
shape of thrust pads. Increasing the aspect ratio of 
elliptical pad (@,) provides higher increment in the 
value of F,, for same strength of magnetic field. 
Percentage increase in the value of F, (based on equa- 
tion (30)) progressively increases from 7.52% for cir- 
cular pad (O) bearing to a maximum value of 14.28% 
for elliptical (E//-IV) pad bearing. 


Nondimensional fluid film stiffness coefficient (S) 


Numerical results for fluid film stiffness coefficient 
(S) of the hydrostatic thrust pad bearing are 


illustrated in Figure 7. These results depict the vari- 
ation of S w.r.t. Hartmann number and geometric 
configuration of pads. Use of electrically conducting 
fluids is noticed to have an adverse effect on the 
stiffness capabilities of the hydrostatic thrust pad 
bearings. The derivate of S w.r.t. H is observed to 
be maximum with the strength of the magnetic field 
lying in the range 3< H <9. In contrast to this, 
the use of the elliptical pads compared to circular 
thrust pads enhances the stiffness capabilities of the 
bearing operating with either of the NL or ECL. 
Hence, changing the bearing pad configuration (i.e. 
use of elliptical pad over circular pad) is observed to 
be compensating the decrease in the fluid film stiff- 
ness occurring owing to the use of electrically con- 
ducting fluids/lubricant. Furthermore, numerical 
result suggests an increase in the value of S as the 
aspect ratio of elliptical pad (œ,) increases. 
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Figure 8. Nondimensional fluid film damping coefficient vs. Hartmann number. 
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The stiffness coefficient obtained for circular (O) 
and elliptical (Ell-IV) hydrostatic thrust pad bearings 
are compared using 


Trendwise: [ 


% S = S lina -5 lentes x woo} (31) 


S | Circular 


H=0—> 10 


The result evaluated after employing the above 
equation indicates that as the strength of the magnetic 
field acting on ECL increases, the percentage incre- 
ment in the value of S also increases. The percentage 
increase in S (based on equation (31)) raised from 
29.77% (H=0) to 44.93% (H=10) owing to the use 
of elliptical (Ell-IV) pad vis-a-vis circular pad in the 
hydrostatic thrust pad bearings. 


Nondimensional fluid film damping coefficient (C) 


Figure 8 shows the variation in the damping coeffi- 
cient of bearing under the combined influence of the 
use of ECL and geometric shape of pads. Results pre- 
sented in Figure 8 suggest that the presence of trans- 
verse magnetic field has a positive and significant 
influence on the fluid film damping coefficient of bear- 
ings. Specifically, this influence is more pronounced in 
strong magnetic field (H > 6). Further, observation 
revealed that invariant to the lubricant employed, 
the use of elliptical pads in the hydrostatic thrust 
pad bearing offers higher fluid film damping as com- 
pared to the circular pads. Increasing the aspect ratio 
of elliptical pad (@,) also enhances the damping cap- 
abilities of the elliptical hydrostatic thrust pad bear- 
ings. Following trends in the curves of S have been 
observed 


[Eleisa = C| enr m Cl ni in = C| we mn < Clena- 


A maximum of 139.37% increase in the fluid film 
damping coefficient has been reported owing to the 
use of elliptical pads and ECL in the hydrostatic 
thrust pad bearings vis-à-vis circular pads. 


Conclusions 


The result and discussion presented in the preceding 
section indicates that the use of electrically conducting 
fluids as lubricant has profound influence on the 
steady-state and dynamic performance characteristics 
of hydrostatic thrust pad bearings. The response of 
the ECL or magnetic field on the performance param- 
eters is observed to be highly dependent on the geo- 
metric shape of the pad. Influence of ECL on the 
steady-state characteristics of hydrostatic thrust pad 
bearings is more pronounced while using elliptical pad 
vis-à-vis circular pad. The performance of elliptical 
hydrostatic thrust pad bearings is observed to be a 
strong function of aspect ratio of elliptical pad (a,). 
Increase in the aspect ratio (skewness) of elliptical pad 
is observed to have an adverse effect on the pocket 
pressure and load-carrying capacity of the bearings. 
Fluid film damping and stiffness coefficient of the 
elliptical pad bearings enhances significantly with an 
increase in the aspect ratio in elliptical pad. Lubricant 
flow requirement (Q) of hydrostatic thrust pad bear- 
ings reduces considerably while operating with elec- 
trically conducting fluid/lubricant in the presence of 
the magnetic field. The circular thrust pad bearing is 
noticed to have lower lubricant flow requirement than 
that of the elliptical thrust pad. The presence of ECL/ 
magnetic field adversely affects the fluid film stiffness 
capabilities of the circular and elliptical hydrostatic 
thrust pad bearings. However, using elliptical shape 
pad in the bearing proved to be beneficial in improv- 
ing the fluid film stiffness coefficient (5) of the bearing. 
On account of using elliptical pad, a maximum of 
44.93% increase in S was observed as compared to 
the circular hydrostatic thrust pad bearings. Fluid film 
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damping properties showed a marked improvement 
due to the presence of ECL and the use of elliptical 
pad in the hydrostatic thrust pad bearing while oper- 
ating with ECL. Increase in the aspect ratio (a@,) of 
elliptical pad further improves the value of damping 
coefficient (C) of the bearing. A maximum of 139.37% 
increase in C has been reported in the elliptical hydro- 
static thrust pad bearings (H = 10; E//-IV) in compari- 
son to the circular hydrostatic thrust pad bearings 
(H=0; O). 
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Appendix 

Notation 

a, b semi-major and semi-minor axis of 

7 elliptical thrust pad (mm) 

a&b aff; b/f 

Ap area of bearing pad (zr?) (mm?) 

Aoc area of circular recess (7 17) (mm?) 

A, Ap/Ave = 4 

B magnetic field vector (Wb/m’) 

B, magnitude of magnetic field vector 
along film thickness direction (Wb/m”) 

C fluid film damping coefficient (N-s/m) 

Č Ch) 

F Pa 

i : (FE) capite 

CS > CS | /6 

E induced electric field vector (V/m) 

Ex, Ez component of electric field along x and 
z directions (V/m) , 

F dynamic fluid film reaction (i Æ 0) (N) 

Fo steady-state fluid film reaction Q = 0) 
(N) 

Fy Lorentz force (N) 

[F] assembled global fluidity matrix 

Fo Fo/reps 

F F/r-ps 

h fluid film thickness (mm) 

ho reference film thickness (mm) 

h h/ho 

H Hartmann number, (H = B,h,(o/ wW! 2) 

J current density (A/m?) 

L capillary length (mm) 

mı, My direction cosines 

M' body couple per unit mass (N-m/kg) 

Ne number of elements 


ny nodes in an element 


tolp 


uU, V, W 


Langrangian interpolation function 
pressure (MPa) 
pocket pressure (%) 4 0 (MPa) 
supply pressure (MPa) 
global nodal pressure vector 
P/Ps 
bearing flow (mm?/s) 
restrictor flow (mm?/s) 
global nodal flow vector 
2 
pn OR 
12 
tabillary tube radius (mm) 
outer boundary radius of circular thrust 
pad (mm) 
radius of the pocket/recess (mm) 
Pofta =1 
r/o = 9.5 
residue of Reynolds equation 
global squeeze matrix 


fluid film stiffness coefficient (N/mm) 
ho 
"Ps 


Time (s) 
Ko 

(i 

tolerance value of convergence for 

nodal pressure vector 

fluid velocity along x-, y-, and z-axes 

(mm/s) 

tangential velocity of runner pad 

(mm/s) 

Cartesian coordinate system 


nondimensional coordinates along x 
and z directions, x/ro; y/fo 

aspect ratio of elliptical pad (major to 
minor axis), a/b 

X/To; Y/To 

material constants for lubricant (N-s) 
viscosity of the lubricant (mPa-s) 
nondimensional angular velocity of 
runner 

MHD functions 

density of lubricant (kg/m*) 
electrical conductivity of lubricant 
(mho/m) 


